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ABSTRACT 

Observations show the presence, in the halo of our Galaxy, of stars moving at velocities 
so high to require an acceleration mechanism involving the presence of a massive 
central black hole. Thus, in the frame of a galaxy hosting a supermassive black hole 
(10® Mq) we investigated a mechanism for the production of high velocity stars, which 
was suggested by the results of iV-body simulations of the close interaction between a 
massive, orbitally decayed, globular cluster and the super massive black hole. The high 
velocity acquired by some stars of the cluster comes from the transfer of gravitational 
binding energy into kinetic energy of the escaping star originally orbiting around the 
cluster. After the close interaction with the massive black hole, stars could reach a 
velocity sufficient to travel in the halo and even overcome the galactic gravitational 
well, while some of them are just stripped from the globular cluster and start orbiting 
on processing loops around the galactic centre. 

Key words: galaxies: haloes - galaxies: nuclei - galaxies: star clusters; stars: kine¬ 
matics and dynamics 


1 INTRODUCTION 

The existence of high velocity stars in the Galactic halo is an 
ascertained feature. Some of them have speed sufficient to 
escape the Galaxy gravitational potential. They may have 
gained such high velocities thanks to different physical mech¬ 
anisms, as three-body interactions among binary systems in 
star clusters or with the massive black hole in the Galactic 
centre. High velocity stars can be divided in two different 
categories, i.e. runaway stars and hypervelocity stars. 

Runaway stars, historically defined in the context of 
O and B stars liHumason fc Zwickvlll94^ . are Galactic halo 
stars with peculiar motions higher than 40 km s“^ (although 
the definition of runaway star is not univocal). Young mas¬ 
sive stars are not expected to be observed in the halo far 
from star-forming regions, since special conditions, as the 
presence of molecular clouds with dense cores, are required 
to form such stars. Therefore, they are thought to be born 
not in the halo, but rather to have travelled far from their 
birthplace. There are two proposed mechanisms for the pro¬ 
duction o f runaway stars: supernova ejections and dynamical 
ejections llSilva fc Napiwotzkill201ll ). _ 

In the supernova ejection mechanism llBlaauwl Il96ll : 
IPortegies Zwartll2000l ) a runaway star is supposed to have 
origin in a binary system when its companion explodes as 
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a supernova. The maximum possible ejection velocity is 
given by the sum of the orbital velocity of the progenitor 
binary and of the supernova kick velocit y. Additional 
effect s may come from asym metric explosions llScheck et al.l 
I2OO6I : IPrzvbilla et al.l l2008l) . but, in any case, runaways 
velocities are below the Galactic e scape velocity. In the 
dyna mical ejection mechanism (IPoveda. Ruiz fc Alien! 

the runaway star derives from a three- or four-body 
interaction. For example, if a binary system interacts with 
a massive star, one member of the binary could be captured 
by the massive sta r, while the other star may be ejected 
with hi gh velocity ([L eonard fc DunccinI 1990|; Gyaramud^ 
2 OO 9 I: iGvaramadze. Gualandris fc Porteeies Zward l2009l : 


Gvaramadze fc Gualandrij 201 iF Perets fc Subr 2OI2I ). In 

this case, the maximum possible ejection velocity is the 
escape velocity of the most massive star. Observations 
show that both the ejec tion mechanisms operate in nature 
(iHooeerwerf et al.ll200ll) . 


Hy pervel o city s tars (HVS) are stars escaping the host 
Galaxy. iHillj (119881) was the first to predict theoretically 
their existence as a consequence of interactions with a 
massive Black Ho le (BH) in the Galactic Centre, while 
iBrown et al.l (l2005l ) serendipitously discovered the first HVS 
in the outer stellar halo, a B-type star moving over twice 
the Galactic escape velocity. The most recent HVS Survey 
is the MMT (Multiple Mirror Telescope) survey, a spectro¬ 
scopic survey of stars within the range of colours of 2.5 — 4 
Mq late B-type stars dBrown et al.l]20l3 ). The MMT ob- 
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servational approach is justified by that such stars should 
not exist at faint magnitudes in the outer halo unless they 
were ejected until those distances. Actually, such stars have 
relatively short lifetimes and should originate in a region of 
on-going star formation. The MMT survey revealed 21 hy¬ 
pervelocity stars ejected from the Milky Way at distances 
between 50 and 120 kpc. However, the MMT survey is able 
to measure only the component along the line of sight of the 
velocity vector. To compute the other velocity components, 
a measure of proper motions is needed. This is i n most of the 
cases not possible (but for some special cases llBrown et al.l 
[2qT3)), because HVSs found by MMT Survey are very dis¬ 
tant so that their proper motions are too small (^1 mas 
yr“^) to be measured with ground-based telescopes. More¬ 
over, as said above, the MMT Survey is biased to the ob¬ 
servation of stars within the colour range of 2.5 — 4 Mq 
late B-type stars, although an observational effort to find 
an older population o f HVS (iBrown. Geller fc Kenvonll2009l : 
iKollmeier et al.l[20091 . l20ldj has been done, unsuccessfully, 
in the last years. 

Hills’ mechanism involves the tidal breakup of a 
binary passing close to a massive BH. This mecha¬ 
nism was analyzed by other authors in the attempt 
to shed light on the propertie s of the stars that 


Gualandris. Porteeies Zwart & Sinioi 20051: Bromlev et al.l 

20061: Sari. Kobavashi & Rossill2009l: iKobavashi et al.ll2012l: 

Rossi. Kobavashi & Sari 20141. Ths 

tidal breakup of a bi- 
on of stars orbiting in 
and the central BH, the 

nary could lead also to a populat 
the inner regions of the Galaxy aro 

so-called S stars (iGould & OuillerJ 

20031: iGinsbure & LoebI 

20061: IPerets. Hooman & Alexander 

20071). Since the Hills’ 


prediction, a lot of other mechanisms have been proposed 
in the literature to explain the production of high velocity 
and hypervelocity stars, which involve different astrophys i- 
cal frameworks and phenomena (iTutukov fc Federoval200^ : 

• the interaction of a Super Massive Black Hole Bi¬ 
nary (SMBHB) with a single star in the nucleus of 
the host galaxy. In this way stars can be accelerated 
at velocities high e nough to escape the local gravita- 


20051: iBaumeardt. Gualandris & Porteeies ZwartI 20061: 

Sesana. Haardt & Madau 20061' 

; 

nother (nearby) galaxy with 
rresent galactic environment 

• some stars can come from a 
a high velocity relative to the r 

IIGualandris & Porteeies ZwartI 

I 2 OO 7 I: iBonanos et al.l 20081: 

Sherwin. Loeb & O’LearvI l2008l: IPeretd l2009l: iBrown et al.l 

2010|); 

• the close encounter of 
star and a single massive 
with velocities larger than 

a hard massive binary 
star could produce stars 
the local escape velocity 


llGvaramadze. Gualandris fc Portegies Zwartll2009l l: 

• a Supernova explosion in a close binary can give a 
kick t o the companion to accelerate it to a very high ve¬ 
locity dPortegies Zwartll2000l : IZubovas. Wvnn fc GualandrisI 
1201311 . 


An interesting point is that , as sh own bv iHanse ^ (l2007li 
and iLopez- Morales fc BonanosI (l2008l l. HVS in our Galaxy 
have both slow and rapid rotations, suggesting different 
acceleration mechanisms. Note indeed that HVSs origi¬ 
nated by a binary disrupted by a Black Hole are not ex¬ 


pected to be fast rotators and so this origin is unlikely 
for fast rotating HVSs. Therefore, since HVS production 
mechanisms should involve different astrophysical frame¬ 
works and phenomena, it would be possible to infer in¬ 
formation about different pieces of physics, as that of the 
three-bod y interaction, the phys i cs of the region near mas- 
sive BHs dGo uld fc QuillerJ [20031 : ISesana. Haardt fc Madaul 
l2007l:lQ’Learv fc Loebll20^’l as well as the physics of Super¬ 
novae. Moreover, the study of the proper motions of such 
fast moving stars can improve the knowledge of the Galaxy 
gravi tational potential shape, of its Dark M atter compo¬ 
nent (iGnedin et al.l[2005l : IYu fc Madaull2007ll and, in line of 

J irinciple, m ay lead to useful information also for cosmology 
Loebfl201l|j . 


Observations of high velocity and hypervelocity ob¬ 
jects have been limited to high-mass, early-type stars 
due to obvious observational bias. Observers have started 
investigating low-mass high velocity star s only recently 
l[PaUadh20_et^jjjj2CI^^ et al.l l2014l : iLi et al.l 120151 : 

IVickers. Smith fc Grebe! 2015l l. some of which are low-mass 
HVS candidates. The European ESA satellite GAIAQ is ex¬ 
pected to measure proper motions with a precision of 0.1 
mas yr“^ and so will be able to provide for a larger and less 
biased sample. Furthermore, Gaia is expected to find ~ 100 
HVSs in a sample of ~ 10® stars. 

The aim of this paper is to investigate another mech¬ 
anism of production of high velocity stars, which involve 
a Globular Gluster (GC) that during its orbit has the 
chance to pass close to an SMBH in the center of its host 
galaxy. This chance is increased by the orbital decay suffered 
by massive clusters moving a dense galactic environment, 
which makes significant the dyna mical friction braking ex¬ 
erted by the stars of the ga l axy (ICapuzzo-Dolcetta 199^; 
iGapuzzo-Dolcetta fc Miocchil l20odr ' Antonini et al.l 2012l l. 
For these test cases, we assumed Mbh = 10® Mq with the 
scope of identify at better the underlying physical mecha¬ 
nism. 

The paper is organized this way: in Sect. [2] we outline 
and describe our approach to the study of the consequences 
of the GG-SMBH interaction; in Sect.|3]the results are pre¬ 
sented and discussed; in Sect. |4] we draw the conclusions. 
Significant details are given in the Appendix. 


2 METHOD 

Our scattering experiments refer to the interaction of three 
different bodies: a super massive black hole (SMBH), a glob¬ 
ular cluster (GG) and a star. In our simulations the SMBH 
sits initially in the origin of the reference frame, while the 
GG follows an elliptical orbit at a relatively close distance 
around it. The assumption of close distance to the BH is 
motivated by that the globular cluster is supposed orbitally 
decayed by dynamical friction braking, as discussed in Area 
Sedda et al. (in preparation). Given the BH influence radius 
as that within which the BH potential dominates 


Tinf — 


GMbh 


( 1 ) 


^ http://www.cosmos.esa.int/web/gaia 
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where a is the stars velocity dispersion in the central 
galactic region, we can treat the dynamics of the star inter¬ 
acting with the environment as a 3-body (GC, SMBH and 
test star around the GC) problem whenever the relevant fly¬ 
by passage occurs within the distance ri^f from the SMBH. 
In the cases studied in this paper (Sect. [3]), < 12.5 pc, 

which means that the choice we make in this paper of 10 pc 
as radius of the GC circular reference orbit, with the values 
of pericenters distances (around which the scattering is ef¬ 
fective) given in Table 1, is fully compatible with the neglect 
of the smooth external field. 

Actually, neglecting the stellar background potential, 
the mechanical energy (per unit mass) of the GC on a cir¬ 
cular orbit of radius is 

1 2 GMbh 1 GMbh 

given that the circular velocity is Vc = [GMbh. 
Consequently, taking into account that the angular momen¬ 
tum per unit mass of the GC on the circular orbit around 
the BH is Lc = \/GMmrH, the pericenter (r_) and apocen- 
ter (r+) distances of the GC on orbits of same energy (Ec) 
but different angular momentum 0 ^ L ^ Lc are given by 



Of course in the above equation the — sign gives the 
pericenter and the -I- sign gives the apocenter. In conclu¬ 
sion, once we have, as reference, a circular orbit of radius 
rc we may compare it with a set of orbits at same energy 
just varying the ratio L/Lc. The eccentricity of the orbit is, 
trivially. 


rjf. — r_ 
r- -t- r+ 



We varied the parameter 


(4) 


a = 



(5) 


in order to sample GC orbits of different eccentricity 
and same orbital energy. Of course, a = 0 for radial orbits 
(e = 1) and a = 1 for circular orbits (e = 0). 

The cartesian reference frame has been chosen as that 
with the a;-axis along the line connecting the GC with the 
SMBH and y-axis orthogonal, so that the (x, y) frame is 
equiverse to the GC orbital revolution. 

In the restricted three-body problem, it is well known 
the existence of t he Hill’s surface s which enclose the two 
hnite-mass bodies (ISzebehelvIl 1966ll . The radius of the Hill’s 
sphere, given by 


XL = ro 


Mac 

3Mbh 


1/3 


( 6 ) 


defines the spherical volume around the GC where its grav¬ 
itational potential dominates. Outside the Hill’s sphere, the 
BH potential overcomes the one of the cluster. For the set 
of parameters used in our scattering experiments, 0.55 pc 


Table 1. The values of a, the eccentricity e, the pericentre r_ 
and apocentre r_|_ of the GC elliptical orbits. 


a 

e 

r_ (pc) 

r+(pc) 

0.1 

0.95 

0.51 

19.6 

0.2 

0.89 

1.06 

18.9 

0.3 

0.84 

1.63 

18.5 

0.4 

0.77 

2.25 

17.8 

0.5 

0.71 

2.93 

17.1 


^ ri, ^ 0.63 pc, 1.19 pc ^ rn ^ 1.36 pc and 2.55 pc 
Si ^ 2.93 pc for a GC of mass 10"^ Mq, 10® Mq and 
10® Mq, respectively. 

A meaningful study refers to the fate of stars moving 
around the GC with orbits initially all within the GC influ¬ 
ence radius. Therefore, we put the initial circular orbits, on 
which the test star moves around the GC, inside this sphere 
by setting the radius of this orbit to be a fraction (< 1) of 
the distance from the first Lagrange point (LI) and the GC. 

For the sake of statistical significance, once fixed the 
unperturbed star circular orbit we sampled cases with initial 
different phases, in the range 0 T 360° at increments of 15° 
llGinsburk Loeb. Wean^l2012l '). 

To summarize, in this paper the values of the relevant 
initial parameters have been set as follows (see also Table 
1 ): 

• the super massive black hole mass is Mbh = 10® Mq; 

• the globular cluster mass, Mac, assumes the three val¬ 
ues 10"*, 10®, and 10® Mq; 

• the test star mass, m*, is set equal to 1 Mq; 

• the GC reference circular orbit has the radius ro = 10 
pc; 

• the GC orbital eccentricity ranges from e = 0.71 (a = 
0.5) to e = 0.95 [a = 0.1) and is parametrized varying 
0.1 ^ a ^ 0.5 at steps of 0.1; 

• the test star circular orbit radius around the GC is 
parametrized by /3 = r/rn, whose values are in the range 
0.08 T 0.25; 

• the star initial position on the circular orbit of given 
radius (see above) is parametrized by adopting a set of 24 
different angles spanning 0 4- 360° with a 15° step; 

• the star circular orbit around the GC and the GC orbit 
respect to the BH are coplanar. 

The choice of the range of a, and consequently of e, 
toward large values of e, is due to that (as we will see in the 
Results Section) the efficiency of the energy transfer on the 
test star orbiting the GC tends to vanish at eccentricities 
less than ^ 0.6. 

Given the above set of initial parameters, we integrated 
the system of the differential equations of the 3-bodies 
(SMBH, GC and star) motion 




,'(ri - rj) 


(7) 


for i = 1,2,3, using the fully regularized algorithm of 
iMikkola fc AarsethI (1200 ih . The need of a regularized algo- 
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rithm is due to the enormous range of variation of the masses 
involved, going from 1 Mq of the test star to 10 ® Mq of the 
SMBH. Any not-regularized direct summation code would 
fail in dealing with the close star-SMBH interaction and 
would carry to an enormous energy error during the close 
triple encounter. 

In the Mikkola’s ARW code this problem is overcome by 
a transformed leapfrog method, which leads to extremely 
accurate integrations of the bodies trajectories when 
combined with the Bul irs ch-Stoer extrapolatio n meth od 
jBuIirsch_&_Stoe3 19^- _ ^kkolu^_&_Thnikawa|_ 1999a 3 


iMikkola fc Merrittll^OOff bOOSl : iHellstrom fc Mikkolall201^ . 


Thanks to the regularized algorithm, the fractional energy 
error is kept below 10“^^ over the whole integration time. 


3 RESULTS 


In our scattering experiments the test star orbiting the GC 
has three possible fates: 

(1) it remains bound to the GC on an orbit significantly 
perturbed respect to the original one; 

(2) it becomes a high velocity star, either bound or un¬ 
bound to the host galaxy; 

(3) it is captured by the massive BH gravitational field 
and starts revolving around it. 


The distinction among these three different situations 
is made by computing the mechanical energy of the star 
respect to the BH and the GC after the scattering. If the 
energy of the star respect to the GC remains negative, the 
star remains bound to the GC, while if this energy becomes 
positive, while the star energy respect to the Black Hole 
is negative, the star becomes bound to the BH. Finally, if 
both these energies are positive, the star is able to leave 
the BH-GC system and, according to the assumed galactic 
potential model, it will be bound or unbound respect to the 
galaxy. The branching ratios, i.e. the probability of different 
outcomes, are plotted in figs. |l(a)|l(b)|l(^ 


Fig. 1(a) gives the branching ratio of stars which are 


captured by the BH after the GC-BH scattering and be¬ 
come bound to the BH, as function of a for the different GC 
masses. As expected, this ratio decreases for larger values of 
a (less eccentric orbits), as well as for larger GC masses. 

Fig. 1(b) shows the branching ratio of stars which re¬ 
main gravitationally bound to the Globular Cluster after 
GC-BH scattering, as function of a for the different GC 
masses. The branching ratio increases for higher values of 
Q, i.e. for less eccentric orbits, and is almost independent of 
the GC mass in the range studied. Note that at values of 
e ~ 0.7 the fraction of bound stars is about 80% for the 10^ 
M© GC and about 67% for the 10® and 10® Mq cases. 

Finally, Fig. 1(c) shows the branching ratio of stars 
which, after GC-BH scattering, leave the BH-GC system 
becoming high velocity stars, as function of a for the differ¬ 
ent GC masses. This fraction decreases for larger values of 
a (less eccentric orbits) and increases for larger GC masses. 

The Lagrangian radius have an important role during 
the close interaction between GC and BH. A star is lost by 
the GC when it crosses the Hill’s surface through or near the 
first or the second Lagrangian point. If it passes through LI, 
its fate is the capture by the BH, while in the second case 


it will escape the whole system, becoming unbound both re¬ 
spect to the BH and the GC. The hrst channel is favoured 
by lower GC to BH mass ratios, since the BH potential is 
stronger and is able to capture a higher number of GC stars 
making them pass through the first Lagrangian point. At the 
same time, the GC gravitational potential is not so intense 
to give the star a velocity high enough to escape the whole 
GC-BH system. Therefore, the branching ratio for the pro¬ 
duction of stars captured by the BH is higher for lower GC 
masses, while the branching ratio of ejected stars increases 
for higher GC masses. 

Another important feature of this mechanism is the sig¬ 


nificant level of collimation of the ejected stars. Figure 2(a) 


shows the average ejection angle (j) as function of a for dif¬ 
ferent GC masses. This angle is that between the velocity 
vector at ejection and the x axis of the inertial reference 
frame, taken in the direction pointing from the BH to the 
initial position of the GC. Figure 2(b) shows the ratio be¬ 
tween the standard deviations and 4>y giving a measure 
of the ejection collimation. If the value of was precisely 
zero, it would mean that all the stars are ejected along the 
same direction after the GC-BH interaction. Therefore, non¬ 
zero values measure the width of the loss-cone within which 
stars are ejected. We see that the average emission angle 
increases for larger values of the GC mass, while the level 
of collimation decreases. Therefore, stars ejected during the 
interaction of a low-mass GC and a massive BH are concen¬ 
trated in a small amplitude jet, which, instead, has a greater 
span for a high ma ss GC. This results are compatible with 
ISesana et al.l (l2006lL who found, studying the HVS produc¬ 
tion of a binary BH, that the collimation is higher when 
the mass ratio between the BHs is lower. Actually, we found 
that stars are highly collimated in jets when Mac = 10^ 
Mq, while for higher GC masses the ejection jets tends to 
become a nearly isotropic emission. 


3.1 High velocity stars 

From our scattering experiments, it is possible to derive the 
velocity profile of the ejected stars. Fig.[3]shows the velocity 
profile of the ejected stars at 20 pc (as it will be justified in 
the following) for different GC masses and all the orbits. The 
distributions are nearly Gaussians, cut at 212 km s“^, which 
is the escape velocity from the BH. Moreover, the peak of 
the distribution depends on the GG mass, the greater the 
velocity peak the larger the GG mass. The ejected stars can 
have two fates: 

(1) to remain gravitationally bound to the galaxy, al¬ 
though having escaped the BH-GC system; 

(2) to become unbound stars, and so HVS, if their kinetic 
energy is sufficient to overcome the galaxy gravitational po¬ 
tential well. 

To evaluate whether stars formerly belonging to the 
GC and ejected at high velocity remain bound to the host 
galaxy, an assumption on the galactic field has to be made. 
In our analysis, we assumed two different models for the host 
Galaxy, one as an elliptical and one as a spiral galaxy. 

The elliptical galaxy potential is represented as a two- 
component model given by a spherical bulge-halo summed 
to the SMBH po tential. The bulge-halo potential is given by 
(lHernauistlll990l ~) 
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0.95 0.89 0.84 0.77 0.71 



a 

Figure 1. Branching ratios of stars captured by the BH (a), GC stars (b) and ejected stars (c), after GC-BH scattering, for Mqq = 10“^ 
Mq (solid line), Mqq = 10^ M© (dashed line), Mqq = 10^ Mq (dot-dashed line) and different GC orbits, parametrized by q = (L/Lc)^. 
The dotted lines represent the branching ratios when the star orbit and GC orbit are perpendicular for Mqc = 10® Mq. 
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Figure 2. Average ‘emission’ angle of the ejected stars (a) and relative dispersion (b), both as function of a, for Mqc = 10^ Mq (solid 
line), Mqc = 10® Mq (dashed line), Mqq = 10® Mq (dot-dashed line). 


^b{r) 


GMb 
r + ab’ 


(8) 


where Mb = 7.8 x 10^° M q and Ub = 5.4 kpc. N ote 
that these values are taken from iMarconi fc Hutiu (l2003l ) to 
represent the elliptical galaxy NGC 3377, whose central BH 
has an estimated mass Mbh = I.OIq ® x 10® Mq. 

Assuming the host galaxy as a spiral, we consider a 
four-component model for its potential 


<E>(r, z) = $fl/f(r) -I- $6(r) -f <f>d(r, 2) -h «E>/i(r), (9) 

where the indexes b, d, h stands for bulge, disk and 
halo, respectively. As above, t he bulge poten tial is ex¬ 
pressed as a Hernquist sphere llHernauistI Il990| ) with the 
constants Mb — 10 ^° Mq a nd q g = 1 kpc, as taken from 
iKornreich fc Lovelace! (l2008l ) and iLingamI (|2014| ). to repro¬ 
duce the paramete rs observed in giant disk galaxies su ch as 
those described by Wang, Sulk anen fc Lovelacd lll992l) and 

iRownd. Dickey &: Heloul 1 19941 ). _ 

_ The axisymmetric disk potential is llMivamoto fc Naeail 

[Hzi) 


^d{r,z) 


GMi 





( 10 ) 


where Md = 10^^ Mq, ad = 6.5 kpc and bd = 0.26 kpc 
llLingamll2014l ). 

Finally, the halo density distribution accounts for the 
presence of a spherical dark matter halo, o f total mass 5.6 x 
10^^ Mq, whose potential is (^Binnev^ll981^ 


^h{r) = ^vlln{r'^+ rl)+c, (11) 

where Vh = 250 kms“^, rn = 2 kpc and c is a constant 
which giv es a match w i th an external (r ^ 50 kpc) keplerian 
potential (lFuiitall2009l : lLingamll2014l ). 


Given a galaxy model, it is possible to calculate the 
escape velocity 


„U) 


{r, z) = ^^-2<Fi(r,2), 


( 12 ) 


which depends, besides on the galaxy model itself 
(j=E,S), also on the position in which it is computed. In 
our 3-body scattering, we follow the trajectories of stars un¬ 
til they are 20 pc far from the BH and compute their velocity 
and the escape velocity at this distance. At this distance the 
contribution of the galaxy components is not negligible in 
the Eq.[T2] but its gravitational potential is nearly constant, 
both for the elliptical and spiral galaxy. Therefore, the mo¬ 
tion of the escaping star is due to only the central BH. 

If u* < Vesc, the star, although escapes the BH-GG sys¬ 
tem, will be gravitationally bound to the galaxy. In this case, 
the star can be bound only to some galactic components, as 
specihed by the vertical lines in Fig.O which divide the dis¬ 
tribution in different portions. The leftmost thick line in all 
the panels indicates the escape velocity from the BH (212 
km s“^). In the left column panels, the other vertical line 
refers to the escape velocity (418 km s“^) from the BH + 
bulge-halo system, while, in the right column panels, the 
other vertical lines indicate the escape velocity respect to 
the BH + bulge (365 km s“^), BH -I- bulge -I- disk (516 
km s“^), BH -I- bulge -I- disk + dark halo (759 km s“^), 
respectively. Therefore, for example, the portion of velocity 
distribution, on the right of the 365 km s“^ vertical line 
(for a spiral galaxy), is unbound respect to the BH -|- bulge 
component but bound respect to the disk and the dark halo. 

On the contrary, if u* > Vesc, it will become a HVS. 
The branching ratios of HVS, i.e. the probability of pro¬ 
ducing unbound stars from the galaxy with respect to the 
total ejected stars, are listed in Tab.[2]for all the orbits. The 
results, which depend on both the total mass of the host 
galaxy and on the shape of its gravitational potential, show 
that the higher is the GG mass the higher is the probability 
of producing HVS. Furthermore, according to the param¬ 
eters chosen for the host galaxies gravitational potential, 
while the 10® Mq GC is able to generate HVS in both the 
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Figure 3. Velocity distribution of escaping stars for Mqq = 10'^ (a,b), Mac = 10^ Mq (c,d) and Mqc = 10^ M© f^ie 

orbits, both for a yitot = 7.81 x 10^*^ Mq elliptical galaxy llMarconi &: Hu^l2003h (left column) and a Mfot = 6.60 x 10^^ Mq spiral 
galaxy llFuiitall2009t) (right column). Vertical lines indicate the escape velocity of the various galactic components (see Sect. [311 


galaxies, the 10^ Mq and 10® Mq GC are able to produce 3.2 The role of star orbital inclination 
HVS only in the elliptical galaxy. 

In our scattering experiments the initial orbit of the star 
and the orbit of the GC are coplanar. In order to check the 
effect of the relative inclination between the star and the 
GC orbit, we performed the same set of simulations for the 




































































8 R. Capuzzo-Dolcetta and G. Fragione 



Figure 4. Comparison between the velocity distributions of es¬ 
caping stars for Mcc = 10® Mq when the GC orbit and star 
orbit are coplanar (solid line) and when they are perpendicular 
(dashed line). The distributions are cut on the left side at 212 km 
s“^, which corresponds to the escape velocity respect to the BH. 

Mac = 10® GC presented above in the case of star orbits 
initially lying on a plane perpendicular to the GC orbital 
plane. The resulting branching ratios are plotted in Fig. 1. 
While the branching ratios of the ejected stars and of the 
stars captured by the BH decrease, the branching ratio of 
stars which remain bound to the GC increases. Therefore, 
the overall effect is that stars tend to remain more bound to 
the cluster in the inclined case respect to the coplanar one. 

Fig. 4 shows the comparison of the velocity profiles. 
The perpendicularity of the orbits makes the distribution 
to peak at lower values of the velocity and the area under 
the distribution is smaller, because the branching ratio of 
ejected stars is lower than in the coplanar case. 

3.3 The role of a smooth GC potential 

In order to see the effect of a GC mass profile in the results, 
we performed the same set of simulations perfor med in the 
case o f a Mac = 10® point mass GC, assuming a lPlummerl 
ImS) mass profile 

M{r) = Mac (13) 

(r2 +a^) I 

where Mgc is the total mass of the GC and a its core 
radius, which is set to 0.5 pc. Fig. 5 shows the velocity pro¬ 
files for a point mass GC and a Plummer GC. The clear 
effect of smoothing the GC potential is that the velocity 
distribution shifts towards lower values of the velocity. Ac¬ 
tually, for the set of parameters chosen in this study the 
gravitational energy of the star is ~ GMcc/a,, which leads 
the peak of the nearly Gaussian distribution to a lower ve¬ 
locity and makes its dispersion decrease. Actually, if the GC 
is taken to be a point mass, for same radius of the circular 
orbit, the generic star of our simulation has a lower (more 
tightly bound) gravitational energy respect to the case of 
a GC smooth potential. Then the amount of gravitational 
energy that could be converted into kinetic energy would be 
higher, giving a larger number of ejected stars and a velocity 



Figure 5. Comparison between the velocity distributions of es¬ 
caping stars for Mgc = 10® Mq when the GC is approximated 
as a point mass (solid line) and when it has a Plummer density 
profile with core radius a = 0.5 pc (dashed line). The distribu¬ 
tions are cut on the left side at 212 km s“^, which corresponds 
to the escape velocity respect to the BH. 


Table 2. Branching ratio (3rd column) of the unbound stars 
(HVS) respect to the total ejected stars for different galaxy mod¬ 
els (E = elliptical, S = spiral). 


Mgc(Mq) 

Galaxy 

BR 

lO'^ 

E 

1.67 X 10-2 

10^ 

S 

0 

10® 

E 

4.47 X 10-2 

10® 

S 

0 

10® 

E 

0.35 

10® 

S 

9.09 X 10-® 


distribution peaked at higher velocities. Clearly, this same 
effect of reduction of the efficiency in the star acceleration 
after the GG-MBH fly-by is obtained when the black hole 
mass is reduced, reducing thus the quantity of gravitational 
energy to inject in the test star motion. This means that we 
would expect a scaling of the phenomenon efficiency almost 
linear with the MBH mass, which would mean a reduction 
of the effects studied in this paper, where the SMBH mass 
is 10® Mq, for a factor 0.04 in the case of the Milky Way 
4x 10® Mq Sgr A* BH. This is just a rough qualitative sketch 
of a context, that of the scattering around the MW central 
MBH, that deserves a much more careful study which we 
will do in a forthcoming paper. 

3.4 Number of ejected stars 

After discussing the fate of an indivudal test star moving 
circularly around a point like GC, we want here to quantify 
the actual number of stars belonging to a GC that can be¬ 
come bound to the SMBH or ejected at high or even hyper 
velocity after the interaction with the SMBH. A suitable es¬ 
timation comes from the evaluation of the number of stars 
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Figure 6. Local fraction G(g) = t'c(C)/r'(€) of stars in nearly cir¬ 
cular orbit with 5 = 0.05, for the Hernquist (a) and the Plummer 
(b) potentials. The isotropic fraction (/3 = 0, solid line) stands 
between the tangentially biased fraction (/3 = —1/2, dashed line) 
and the radially biased fraction (/3 = +1/2, dot-dashed line). 


in nearly circular orbits in a self consistent model of GC of 
known distribution function (DF). Of course, the DF is not 
unique and depends both on the functional form of the GC 
gravitational potential and on the level of anisotropy. Here 
we consider two different gravitational potentials for the GC 

• a Hernquist potential (lHernauistlll990l ) 


where (Jr,oe, 0 - 4 , are the velocity dispersions in spherical 
polar coordinates, an isotropic model has /3 = 0, while /3 is 
non-zero for anisotropic cases. 

Using different modelisations, we computed the local 
fraction of stars of the GC on nearly circular orbits, G(^) = 
Vc(H)/vis,), as function of ^ = r/a, where u{^) is the local 
density of bound stars of any velocity, while is that of 

stars moving on nearly circular orbit of scaled radius i.e. 
those with velocity very close to the local circular velocity 

Vc(^). 

For the Hernquist potential. 




M 1 

2iTa^ ^(1 + ’ 


while for the Plummer potential 


(17) 


47ra3 (l + C=)®/2- 

The function vdO depends, besides the assumed po¬ 
tential, also on the degree of velocity anisotropy and can 
be calculated as described in Appendix. The resulting i'c(C) 
will depend also on the ‘tolerance’ S, which quantifies the 
departure from the exact circular velocity. In our calcula¬ 
tions we consider in ^he stars having a local speed 

in the interval —J ^ v/vc ^ +<5, with 5 = 0.05. 

As expected. Figs. |6(a)|6(b')] show that the isotropic 
(/3 = 0) model fraction stays between the radially biased 
{P = +1/2) and the tangentially biased {P = —1/2) models. 

Gombining these evaluations of the local fractional 
abundance of stars on nearly circular orbits with the branch¬ 
ing ratios obtained in the previous sections, we can evaluate 
the actual numbers of stars which escape from or remain 
bound to the GG or the SMBH after close GG-SMBH in¬ 
teractions. This requires some assumptions. First of all we 
assume that the GG is composed by a single-mass (m* = 1 
Mq) population of stars. Moreover, we assume that the core 
radius has a size which is half the innermost circular orbit 
radius chosen. 

The number of stars in nearly circular orbits is evalu¬ 
ated by mean of the integral 


Nc = iirAf 



(19) 


$(r) 


GM 
r + a ’ 


(14) 


where M is the total G G mass and a it s core radius; 
• a Plummer potential (IPlumme 3[T9i3) 


$(r) 


GM 
y/r^ + 


(15) 


where, again, M is the total GC mass and a its core radius. 


Subsequently, we treated both the isotropic and 
anisotropic cases. With the usual definition of the anisotropy 
parameter 


P = l- 


„2 I „2 

2 a^ 


(16) 


where Af = 10^, 10®, 10® is the assumed number of stars 
in the GCs, while Va and ri, are given by 


' a,b 


GM{r 


( 20 ) 


being Va,b = Uc ± Svc- A straightforward product of Nc 
with the fractions of bound to the GC or BH or unbound 
(high velocity or hypervelocity) stars, as derived by our scat¬ 
tering experiments presented in the previous sections, gives 
an estimate of the actual number of stars ejected thanks to 
the GC-BH interaction mechanism. 

Figures 7-8-9 report the number of stars ejected from 
the GC after the close interaction between the GC and the 
BH, as function of the circular orbit radius for various values 
of a, both in the isotropic (P = 0) and anisotropic (P = 
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Figure 7. Number of stars ejected by the Globular Cluster with mass Mqc = lO'^ Mq as function of for Hernquist (left) and Plummer 
potential (right), for different values of the anisotropy parameter 0. The different line styles refer to different GC orbit: a = 0.1 solid 
line, a = 0.2 dashed line, a = 0.3 dot-dashed line, a = 0.4 and a = 0.5 orbits give no ejected stars. 


±1/2) Hernquist and Plummer GC models. The number of 
ejected stars depends, besides on the GC mass, on the GC 
potential model and on the degree of anisotropy, as shown 
by Fig. 6. 


Moreover, the number of ejected stars increases with the 
GC mass. Actually this mechanism would produce 10^ 10^, 
10^ ^ 10^ and 10^ -b 10^ high velocity stars for GC mass of 
lO'^ Mq, 10® Mq, 10® Mq, respectively. 
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Figure 8. Number of stars ejected by the Globular Cluster with mass Mqc = 10^ Mq as function of for Hernquist (left) and Plummer 
potential (right), for different values of the anisotropy parameter 0. The different line styles refer to different GC orbit: a = 0.1 solid 
line, a = 0.2 dashed line, a = 0.3 dot-dashed line, a = 0.4 dotted line, a = 0.5 double dot-dashed line. 


This almost linear dependence of the number of ejected 
stars on the GC mass is explained by that the gravitational 
energy per unit mass available to be converted into the test 
star kinetic energy is linearly scaling with Mac, Egr ~ 
Mac- 


The same token holds for an explanatioir of the 
isotropic-anisotropic models difference. In the case of tan¬ 
gentially biased models (/3 = —1/2) there is a larger frac¬ 
tion of ejected stars because of the larger fraction of stars 
in irearly circular orbits. On the other hand, the fraction 
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Figure 9. Number of stars ejected by the Globular Cluster with mass Mqc = 10® Mq as function of for Hernquist (left) and Plummer 
potential (right), for different values of the anisotropy parameter (3. The different line styles refer to different GC orbit: a = 0.1 solid 
line, a = 0.2 dashed line, a = 0.3 dot-dashed line, a. = 0.4 dotted line, a. = 0.5 double dot-dashed line. 


of ejected stars is smaller for radially biased models (/3 = 4 CONCLUSIONS 

+ 1 / 2 ). 

The phenomenon of the existence of high velocity stars, or 
even hyperveloci ty stars, in our Galax y has been explained 
in the literature (IYu fc Tremainell2003l l in terms of the pres¬ 
ence of a massive black hole in the galactic centre (~ 4 x 10® 
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Mq in our Milky Way). The mechanism of star acceleration 
to high speed requires, indeed, an efficient energy exchange 
in a multiple system, i.e. a t hree or mo re bodies interaction 
with the massive black hole (lHillg|l98f ) or, if it exists, with 


a black hole binary dSesana et al.l 


20091 . 


In this paper we deepened what has been recently found 
by Area S edda et al. (in prep aration) and preliminarily pre¬ 
sented in (ISnera et al.ll2015l l. i.e. that the close passage of 
a massive globular cluster near a massive black hole can be 
source of ejection of stars from the cluster, which are accel¬ 
erated to high speeds. In our study, we assumed Mbh ~ 10® 
Mq with the scope of identify at better the actual physical 
mechanism. The underlying mechanism is likely a 3-body 
interaction, where the ‘bodies’ are the super massive black 
hole (10® Mq), the globular cluster (10^, 10® and 10® Mq) 
and the test star (1 Mq) belonging to the globular cluster. 

We have performed a series of high precision integration 
of the orbits of such 3 bodies to check the probability for 
the test star orbiting a globular cluster, which experiences 
a close (pericenter distance less than 10 pc) encounter, to 
remain bound to the cluster, to be captured by the black 
hole or gain high velocity such to overcome the cluster es¬ 
cape velocity and, possibly, the galaxy escape velocity. We 
determined the branching ratios of these 3 phenomena and 
found that: 


• the efficiency of the star acceleration process is almost 
linear in Mac', 

• a massive globular cluster (composed by 10® identical 
1 Mq stars) releases, in a single close passage around the 
super massive black hole, about 10^ stars; 

• in a very close GC-BH encounter [Mac = 10® Mq, 
a = 0.1) the probability of stars to remain bound, become 
bound to the BH or escape from the cluster are ~ 5%, ~ 
45%, ~ 50%, respectively; 

• the fractions of stars, respect to the total ejected stars, 
which escape from the whole galaxy is ~ 18% for an Mtot = 
7.81 X 10^® Mq elliptical and ~ 0.5% for an Mtot = 6.60 x 
10^^ Mq spiral galaxy. 

Moreover, we studied the effects of the inclination of 
the star initial orbit around the GC respect to the GG or¬ 
bit. When the orbits are perpendicular, the branching ratios 
of the ejected stars and of the stars captured by the BH 
decrease, while the branching ratio of stars which remain 
bound to the GG increases. Therefore, the overall effect of 
the increasing inclination of the star orbit respect to the 
GG orbit consists mainly in a slight reduction of the frac¬ 
tion of stars ejected from the GG and/or captured by the 
BH. Correspondingly, the velocity profile of the ejected stars 
is peaked at lower velcoity values. 

Furthermore, we performed the same set of 3-body scat¬ 
tering experiments for the Mac = 10® GC assuming it 
is a Plummer sphere with a core radius of 0.5 pc. In this 
smoothed case, the velocity distribution shrinks towards 
lower values of the velocity since the amount of gravita¬ 
tional energy, to convert into kinetic one, decreases. As a 
consequence, we found no ejection of hypervelocity stars, 
suggesting that, when the GC is described by a smoothed 
potential, it is likely to eject only high velocity stars, which 
escape the GC-BH system, but still remain bound to the 
host galaxy. 

Different high velocity and hypervelocity ejection mech¬ 


anisms predict different spatial distributions, velocity dis¬ 
tributions an d physical an d kinematic characteristics of the 
ejected stars (lBrownll2015l 'l. 

One important feature of this mechanism is the colli- 
mation of the ejected stars. We found that the high velocity 
stars are ejected in sort of jets, whose angular amplitude de¬ 
pends on the GC mass. High velocity stars produced by the 
interaction of a low-mass GC with a massive BH are likely 
to be concentrated in a small amplitude jet. On the other 
hand, higher GC masses make the ejection jets have a huge 
amplitude and stars are ejected in a n early isotropic emis - 
sion. These results are compatible with lSesana et akl (l2006l l 
finding. Note that our results refer to a single GC-BH inter¬ 
action; anyway, we expect a production of a high velocity 
star jet every time a GC undergoes a close encounter with 
a massive BH. Consequently, small clusters of high velocity 
stars are expected to be present in the sky, whose charac¬ 
teristics depend strongly on the mass of the GC from which 
they come from. At this regard, it is interesting to note that 
about half of the discovered HVS ar e clumped around the 
Leo constellation llBrown et al.|[20l3 i. The kinematics and 
dynamics of the jets, their position in the sky and their 
velocity profile may give information about the GC-BH in¬ 
teraction that produced them, in particular about the GC 
mass and its orbits before being disrupted by the massive 
black hole. Moreover, some stars of the jet, the ones on the 
right tail of the velocity distribution, may have such high 
velocities to be HVS, and so to be lost by the jet itself since 
they are unbound respect to the hosting galaxy gravitational 
field. Being composed by stars formerly belonging to a GC, 
we expect that the jets contains stars of about same age and 
metallicity, with a common flight time. Hence, in principle, 
by measuring the dynamical and physical properties of such 
high velocity jets, it would be possible to infer information 
about the globular cluster progenitor before its tidal erosion 
by the massive BH. 

Finally, a counterpart of such jets is the production of 
a population of stars orbiting the inner galactic regions and 
which share age and metallicity with the jet stars. 

As mentioned above, we assumed Mbh = 10® Mq with 
the scope of identify at better the underlying physical mech¬ 
anism. For what concerns the Milky Way, we expect that 
the velocity distribution would be peaked at lower velocities 
due to the less energetic mechanism (being the SgrA* mass a 
factor 25 lower), with the production of high velocity rather 
than hypervelocity stars. However, we suggest that, in the 
very last, eccentric and narrow GG orbits, some HVS may be 
produced through the mechanism studied in this work. Fi¬ 
nally, we suggest that the recent observe d runaway RR Lyrae 
variable star, MAGHO 176.18833.411, dKunder et akllioisl l 
may have been produced through the mechanism studied in 
this work. 
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APPENDIX A: FRACTION OF PARTICLES IN 
NEARLY CIRCULAR ORBITS 

A spherical stellar system confined by a steady spherical 
potential <I?(r) can be described by a unique ergodic dis¬ 
tribution function (DF). The DF depends on the phase- 
space coordi nates only through the Ha miltonian R(r, v) = 

+ 'I?(r) dBinnev &: Tremaine! l201lh and can be written 
as a non negative f{S), where f = — R J- 4?o is the relative 
energy, with a constant chosen such that / > 0 for £ > 0 
and / = 0 for £1 0. 

The functional expression of f{£), which obviously de¬ 
pends on the functional form of 4?(r), can be found by the 
Eddington’s inversion formula 
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where v(r) is the spatial number density generated by 
the DF 


rf^v/(r, v) = 47ry dv v'^- -n^), (A2) 

and where = —<I> + $o- _ 

The Hernquist potential llHernauistI Il990l i, for which 
"ho = ‘F(r —>• oo) = 0 , 
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The expression of fi{S) depends on the form of the 
potential $(r) and on the value of the anisotropy parameter 
/3. In the case fd = +1/2 the expression of fi{£) is given by 


MS) 

while in the case 


1 _d_ 
2 'k‘^ dT 
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where M is the total mass and a is the core radius, 
leads, by Eq. lAll to 


f{£) = 


Vi 


v/2(27r)3(GMa)(3/2) (1 - T)2 

3 arcsin x/F 


(l-2f)(8f -8f-3) + 
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where £ = £J (GMJa) is t he adimensional relative en¬ 
ergy. The Plummer potential (IPlumme 3 I 19 IIII , for which 
$0 = ‘I>(r —>• 00 ) = 0, 


$(r) 


GM 

a?' 


(AS) 


where, again, M is the total mass and a is the core 
radius, leads, by Eq. lAll to 


f{S) = 


96g'^/^ 

7V8Tv^(GMaf3/2)’ 


(A6) 


where £ is, again, the adimensional relative energy. 
The above DFs are isotropic, i.e. the anisotropy param¬ 
eter 


P = l- 


2 I 2 

2 a^ 


(A7) 


is null. In Ea. lA7l ar, ae, are the velocity dispersions 
on the r, 0 and (j> component, respectively, in a spherical 
polar reference frame. 

Models with /3 yf 0 can be generated by taking the 
distribution function in the form dBinnev fc Tremainell201ll l 


f{£,L)=L-^^h{£). (A8) 

In Ea. lA8l L is the absolute value of the specific angular 
momentum and fi(£) is an arbitrary non-negative function 
of £. Given the DF in the form of Eg. lASI the space number 
density is given by 
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In the case of the Hernquist potential (IA3I) we thus have 
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for P = +1/2, and 
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for P = -1/2. 

The Plummer potential instead, leads to 
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A£^ - ® 


Sn^GMb (1 _ £2^/2 ’ 
for P = +1/2, and to 


(A14) 


MS) = 


30T'‘ - 47T® + 20T® 


8 tt^ (GMb)^ 


(1 _£:2)l/2 


(A15) 


for P = -1/2. 

Finally, to calculate the fraction of stars on nearly cir¬ 
cular orbits G{p) = Vcpp)/v(:p), we consider a ‘tolerance’ <5, 
which quantifies the departure from the exact circular ve¬ 
locity. In our calculations we consider in VciP) all the stars 
having a local speed in the interval —5 ^ v/vc 4 +d. Con- 
seqeuently, their number density is given by 


/■ 


j/c(r) = ldvf(r,v) = 

rvc+Svc 


= h 


where 


f 

«/ Vr- 


dv w^/(4' — 


(A16) 


Ib = 


47r 


p = o 


2 n Jf drj sin 77 P = ± 1 / 2 . 


(A17) 


This paper has been typeset from a T|/X/ DTj^X file prepared 
by the author. 






























